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INTRODUCTION 


There  is  no  optical  region  more  amenable  to  theoretical  treatment 
of  nonspherical  and  compound  particles  than  the  Rayleigh  region.  Here 
it  is  possible  to  calculate  with  relative  ease,  the  extinction  due  to 
any  particle  which  can  be  represented  by  an  ellipsoid^1)  or  a  con- 
focally  coated  ellipsoid^),  This  is  not  very  restrictive  because 
nearly  all  convex  particles  can  be  accurately  described  by  an  ellip¬ 
soid  envelope,  with  the  exception  of  polyhedra.  Coated  aphereB,  need¬ 
les  and  discs  are  accurately  represented.  Recently  lattice  dynamical 
and  dielectric  continuum  calculations  were  made  predicting  extinction 
by  Rayleigh  cubes  and  polyhedra'3“'3' .  It  is  even  possible  to  model 
accurately  the  extinction  caused  by  a  collection  or  irregular  shapes, 
the  most  common  form  for  solid  aerosols.  From  the  standpoint  of  pro¬ 
ducing  strong  extinction  per  unit  mass,  the  Rayleigh  region  is  the  most 
fruitful  for  even  a  moderately  absorbing  particle.  Here  extinction  as 
a  function  of  size  attains  a  high  plateau  which  remains  independent  of 
size  so  long  as  we  stay  within  the  region.  Here  we  do  not  have  to 
depend,  bb  is  often  the  case  elsewhere,  on  narrow  extinction  resonances 
which  all  but  disappear  when  particle  size,  shape,  composition  or 
orientation  change  only  slightly.  The  Rayleigh  region  holds  the  most 
promise  for  yielding  particles  engineered  to  produce  strong  broad 
electromagnetic  extinction,  and  therefore  it  should  be  explored  as 
extensively  as  possible. 

DISCUSSION 

The  Rayleigh  region  is  constrained  by  definition  to  be  located 
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where  particle  size  is  small  wibh  respect  to  wavelength  both  inside 
and  outside  the  particle.  Here  there  is  a  close  connection  linking 
the  static  polarization  of  a  particle  and  its  extinction  per  unit  mass. 
Before  describing  this  connection  first  we  will  take  a  closer  look  at 
the  extinction  per  unit  mass  which  will  be  referred  to  simply  as  the 
extinction.  Extinction,  determines  the  fraction  of  incident 
radiation,  T,  which  passes  a  distance,  L,  through  an  aerosol  cloud 
having  concentration,  C,  in  accordance  with  Beers  Law. 
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A  convenient  self-consistent  Bet  of  units  puts  extinction  in  square 
meters  per  gram  of  aerosol,  concentration  in  grams  of  aerosol  per 
cubic  meter  of  air,  and  pathlength  in  meters.  The  extinction  co¬ 
efficient  depends  on  the  geometric  cross  section,  G,  optical  extinction 
efficiency  factor,  Q,  and  particle  weight,  W,  in  the  following  way. 
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The  inner  brackets  represent  an  average  over  solid  angle  to  take  into 
account  random  orientations  experienced  by  particles  in  the  cloud 
while  the  outer  brackets  represent  an  average  over  particle  mass  (size) 
distribution  in  the  cloud.  This  double  integral  simplifies  to  the 
following  single  integral  expression  because  <GQ/W>  is  independent  of 
particle  size  for  an  absorbing  particle  in  the  Rayleigh  region  and  of 
course  particle  weight  is  independent  of  orientation. 
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The  fundamental  extinction  theorem  relates  the  extinction  cross  section 
GQ,  to  the  real  part  of  the  scatter  amplitude  in  the  forward  direction, 
Rejs(O)} ,  for  radiation  of  wavelength, 


GQ  -  “Re{s(0)} 


The  Rayleigh  theory  relates  the  static  complex  polarizability,  a ,  of  an 
absorbing  particle  to  its  forward  scatter  amplitude. 

S(0)  -  i  ^~a 

Substituting  this  value  into  the  fundamental  extinction  theorem  we 

find  for  a  single  particle  at  one  orientation, 

2 

Gq  .  i-I  Re  {(la )} 

The  extinction  cross  section  resulting  from  a  collection  of  randomly 
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oriented  ellipsoidal  particles  oan  be  expressed  as  a  function  of  a 
complex  polarizability  which  is  just  one  third  the  sum  of  the  polari¬ 
zabilities,  a j,  along  the  major  ellipsoidal  axes. 

a  “  (a  +  a  -tot  )/3 

1  2  3 

The  major  axle  polarizabilities  of  an  ellipsoidal  particle  are  directly 
proportional  to  particle  volume,  V,  and  depend  on  complex  refractive 
index,  n,  and  depolarization  factors,  L j ,  along  the  major  ellipsoidal 
axes  a,  b,  and  c.  , 

V (n  -1) 

“j  “  Lj  (nz  -  1)  4  1 

where  n  ■  N  +  iK 

/“  abc  ds 

2  TS+Ppn  (i+bTjm  (,+c!Z)i n 

with  appropriate  cyclical  changes  for  L2  and  L3.  The  three  depolari¬ 
zation  factors  have  the  properties  that  each  is  nonnegative  and  their 
sum  is  one. 


The  depolarization  factors  of  spheroids  depend  only  on  the  aspect 
or  axial  ratio  of  the  axis  of  symmetry  with  respect  to  either  one  of 
the  other  two  equal  major  axes. For  a  prolate  spheroid  with  axial  ratio, 
A,  defined  always  to  be  greater  than  one; 
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l-2; 


L2  -  1-  Li  L3  -  L2  e2  -  1-1/A2 

For  an  oblate  spheroid  with  aspect  ratio,  A,  also  defined  to  be  greater 
than  one ; 

1l“  1 +--f  Cl  -  I  arc  tan  f) 
f  2  f 

l2"  1  ~  Li  L3  -  L2  f2  -  A2  -  1 


We  may  now  write  the  extinction  for  a  cloud  of  randomly  oriented  Ray¬ 
leigh  ellipsoidal  or  spheroidal  particles. 


^  m  8tt  3  V  2NK 

~3XW  Zj„1  “ {  Lj  (X2-k2)"  +  (1-Lj  ) p  +  ,( 2NKLj ) 
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When  the  wavelength  and  density,  p,  are  moved  over  to  the  left  hand 
side  of  the  equation,  shape  dependence  is  more  clearly  evident  on 
the  right  hand  side.  In  this  way,  wavelength  and  density,  which  are 
altogether  unrelated  to  shape,  no  longer  appear  as  independent  variables 
and  extinction  values  will  not  be  tied  to  a  single  wavelength  and 


density. 
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3  j-1  (NZ-K*)  +  (1-Lj)|T+  (2NKLj)' 


it  is  this  quantity ^pX  that  we  choose  to  plot  on  the  contour  maps. 

This  frees  us  to  select  any  wavelength,  density  combination  and  inter¬ 
pret  extinction  isopleths  accordingly .  One  self-consistent  set  of 
units  convenient  to  uBe  in  the  infrared  holds  the  extinction  in  square 
meters  per  gram,  puts  density  in  grams  per  cubic  centimeter,  and  wave¬ 
length  in  micrometers.  In  this  system  of  units,  extinction  density 
wavelength  contour  maps  may  be  said  to  represent,  for  example,  extinc¬ 
tion  at  one  micrometer  and  unit  density  or  at  half  a  micrometer  and  a 
density  of  two. 


When  the  shape  and  dielectric  properties  of  the  particle  combine 
to  reduce  the  value  of  the  denominator  in  the  previous  equation,  there 
is  an  increase  In  extinction  or  equivalently  a  resonance ,  The 
shape  dependence  of  this  resonance  enters  through  the  depolarization 
factor.  A  sphere  has  depolarization  factors  equal  to  1/3  for  fieldB 
applied  along  any  orthogonal  set  of  three  radial  axes.  A  prolate 
spheroid  approximating  a  thin  needle  has  L^Q  for  fields  applied 
parallel  to  its  length  and  Lj-*'  1/2  for  fields  in  the  plane  of  symmetry. 
An  oblate  spheroid  approximating  a  thin  disk  has  Lj+1  for  fields 
applied  perpendicular  to  the  plane  of  symmetry  andJLj-*0  for  fields 
in  the  symmetry  plane. 

A  resonance  occurs  when  the  following  conditions  are  satisfied! 

Lj  (N2  -  K2)  +  ( 1-L j )  -  0  2NKLj  -  0 

Remembering  that  all  depolarization  factors  are  nonnegative,  we  find 
that  for  a  given  shape  and  depolarization  factor  the  resonant  values 
for  the  optical  constants  which  satisfy  both  conditions  are 

N  -  0  K  -  (1/Lj  -  1)1/2 

A  close  look  at  extinction  in  the  limit  n-M3  and  K  equal  to  the  above 
value  reveals  that  this  in  indeed  an  extinction  pole.  All  three 
depolarization  factors  for  a  sphere  are  equal  to  1/3  and  its  resonance 
therefore  occurs  at  N*0  and  K“  ^2. 
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The  special  case  where  Lj**0  warrents  some  discussion.  This  can 
occur,  as  mentioned  before,  for  selected  orientations  of  needles  and 
disks  with  respect  to  the  electric  field  vector.  Taking  this  limit 
in  the  previous  equation  relating  extinction  to  complex  refractive 
index  and  depolarization  factor,  the  terms  in  the  summation  corres¬ 
ponding  to  depolarization  approaching  zero  become  simply 

tfpA  .  l&L  £  NK  +  terms  with  L.j*  0 

J  Lj-0  J 

These  terms  tend  to  dominate  over  terms  which  correspond  to  depolari¬ 
zation  which  is  nonzero  over  most  of  the  complex  refractive  index 
plane. 

Figure  1  shows  extinction  density  wavelength  contours  for  a  sphere 
over  a  complex  refractive  index  range  typical  of  metals  at  visible 
wavelengths  and  semiconductors  at  visible  and  infrared  wavelengths. 
Future  references  to  extinction  density  wavelength  contours  will  be 
made  simply  to  extinction  for  brevity.  There  is  one  resonance  or 
region  of  strong  extinction  centered  at  N*0  and  K" /2.  This  has  been 
Identified  as  the  first  electrostatic  surface  polariton  mode.  The 
region  where  the  imaginary  component  of  the  refractive  index  is  greater 
than  the  real  component  is  the  restrahl  region.  There  must  always  be 
absorption  here  in  order  to  satisfy  the  KramerB-Kronig  relationships. 
Any  point  lying  on  the  K  axis  violates  this  requirement.  Therefore 
the  K  axis  may  only  be  approached,  and  extinction  will  always  remain 
bounded.  The  closest  approach  to  this  resonance  will  be  made  by  a 
material  with  maximal  oscillator  strength  in  a  single  oscillator 
within  sum  rule  limits  and  with  minimal  oscillator  damping^?). 

In  figure  2  the  Bingle  resonance  of  a  sphere  has  split  into  two 
resonances  for  an  oblate  spheroid  at  an  aspect  ratio  6£  ten.  One  pole 
moves  down  the  K  axis  toward  the  origin  without  noticeable  changes  in 
the  extent  of  its  contour  values  while  the  second  resonance  climbs 
steadily  up  the  K  axis,  spreading  its  influence  in  the  form  of  high 
extinction  over  more  and  more  of  the  complex  index  plane  while  aspect 
ratio  grows.  A  similar  situation  applies  to  prolate  spheroids  except 
that  the  resonance  moving  down  the  K  axis  converges  on  the  value  one. 

A  metal  is  the  aerosol  material  best  suited  to  take  advantage  of 
high  extinction  produced  by  the  resonance  moving  up  the  K  axis.  A 
simple  model''  for  the  optical  constants  of  a  metal,  the  Drude  theory 
model,  puts  the  real  parts  of  the  refractive  index  equal  to  the  imag¬ 
inary  index  at  wavelengths  greater  than  ten  microns.  Extinction  pro¬ 
duced  by  metal  prolate  spheroids  obeying  the  Drude  theory  appears  in 
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Figure  1  Extinction  lsopleths  for 
a  Rayleigh  sphere  as  a  function  of 
the  complex  refractive  Index. 


Figure  3  Extinction  lsopleths  for 
a  metal  (N-K)  Rayleigh  prolate 
spheroid  as  a  function  of  the  base 
ten  log  of  the  aspect  ratio  and  the 
imaginary  part  of  the  complex  re¬ 
fractive  index. 


Figure  2  Extinction  lsopleths  for 
a  Rayleigh  oblate  spheroid  with  an 
aspect  ratio  of  ten  as  a  function 
of  complex  refractive  index. 


Figure  4  Extinction  difference 
lsopleths  between  metal  (N"K) 
spheroids,  prolate  minus  oblate 
compared  at  equal  aspect  ratio 
and  complex  refractive  index. 
Negative  isoplethB  are  deleted. 
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figure  3  as  a  function  of  the  base  ten  logarithm  of  the  aspect  ratio 
and  optical  constants.  This  range  of  optical  constants  Is  typical  of 
metals  in  the  infrared.  Metal  oblate  spheroids  generate  a  similar 
appearing  set  of  contours.  In  figure  4  extinction  produced  by  a 
metal  oblate  spheroid  is  subtracted  from  that  produt eu  d”  ■■  prolate 
spheroid  having  the  same  aspect  ratio  and  optical  c'’-'stants .  Negative 
contour  values  are  not  plotted.  In  figure  5  prolate,  ip ceroid  extinc¬ 
tion  is  subtracted  from  oblate,  so  the  region  where  an  oblate  spheroid 
is  superior  to  a  prolate  spheroid  appears  in  the  area  which  is  con* 
toured  out  to  an  aspect  ratio  of  ten  thousand.  Increasing  metallic 
refractive  index  and  increasing  aspect  ratio  work  to  produce  high 
extinction  as  can  be  seen  in  figure  3.  When  comparing  which  shape 
causes  higher  extinction  at  a  given  complex  index  we  see  from 
figure  4  and  3  that  the  prolate  shape  1b  superior  for  aspect  ratios 
less  than  some  value  while  the  oblate  1b  superior  for  aspect  ratios 
greater  then  that  value.  In  comparisons  at  a  fixed  aspect  ratio,  we 
see  that  the  oblate  spheroid  in  superior  below  certain  complex  indices 
and  for  higher  complex  indices  a  prolate  spheroid  is  better. 

The  original  ellipsoidal  theory  of  Rayleigh  has  been  extended  in 
recent  years  to  include  confocally  coated  ellipsoidal  particles.  The 
major  axis  polarizabilities,  at  ,  of  a  confoca.lly  coated  ellipsoid 
are  directly  proportional  to  overall  ellipsoid  volume,  V,  and  depend 
on  inner  (core)  complex  refractive  index,  ni,  outer  (shell)  complex 
refractive  index,  n0,  core  volume  fraction,  Vi/V,  inner  depolarization 
factor,  L^,  and  outer  depolarization  factor,  Lq. 

..  .JL  +  v}+  CLJ  u'no!)+  Vj  'i 

Where  the  depolarization  factors  and  complex  indices  are  defined  as 
before  for  solid  ellipsoids  and  spheroids.  The  confocal  constraint 
requires  that  the  difference  between  the  square  of  the  outer  axes  and 
the  square  of  the  inner  colinear  axes  remain  constant,  resulting  in 
a  thinner  coating  on  larger  axes.  The  value  for  this  constant  is 
determined  by  the  core  ellipsoid  volume  fraction.  It  is  important  to 
recognize  that  a  confocal  constraint  removes  a  degree  of  freedom  from 
simultaneous  specification  of  inner  depolarization  factor,  outer 
depolarization  factor  and  core  volume  fraction.  Once  any  two  of  these 
variables  are  specified,  the  third  is  automatically  determined  by  the 
constraint. 

If  the  coated  ellipsoid  iB  a  coated  prolate  spheroid,  both  the 
inner  and  outer  spheroids  will  be  prolate  with  two  foci  pointB  In 
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Figure  5.  Extinction  difference  isopleths  between  metal  (N*K) 
spheroids!  oblate  minus  prolate,  compared  at  equal  aspect  ratio  and 
equal  complex  refractive  index.  Negative  IsoplethB  are  deleted. 


Figure  6.  Extinction  difference  isopleths  coated  minus  homogeneous, 
between  a  metal  mantle  (N*20,  K«*20)  dielectric  core  (N*l,5,  K».01) 
oblate  spheroid  and  a  homogeneous  oblate  spheroid  of  the  same  metal. 
Aspect  ratios  of  the  core  and  homogeneous  spheroid  are  equal. 
Negative  contours  are  deleted. 
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Figure  7.  Extinction  difference  isopletba,  coated  minus  homogeneous, 
between  a  metal  mantle  (N-20,  K-20) ,  dielectric  core  (N-1.5,  K-.01) 
oblate  spheroid  and  a  homogeneous  oblate  spheroid  of  the  same  metal. 
Aspect  ratios  of  the  mantle  outer  surface  and  homogeneous  spheroid 
are  equal.  Negative  contours  are  deleted. 


Figure  8.  Extinction  isopleths  for  a  metal  mantle  (N-20,  K-20)  oblate 
Bpheroid  at  an  outer  aspect  ratio  of  100  and  a  core  volume  fraction  of 
0.5  as  a  function  of  core  complex  refractive  index. 


I  .£,-»*  -•*  .’Wi 


EMBURY 


common,  and  if  It  Is  a  coated  oblate  spheroid  then  both  will  be  ob¬ 
late  with  a  foci  ring  of  points  in  common.  Both  inner  and  outer 
depolarization  factors  depend  upon  aspect  ratios  as  they  did  for  solid 
spheroids.  However,  due  to  a  constraint  a  choice  must  be  made  as  to 
which  two  of  the  variables,  Inner  aspect  ratio,  outer  aspect  ratio  or 
core  volume  fraction  will  be  specified.  Inner  aspect  ratio  is  always 
greater  than  or  equal  to  outer  aspect  ratio.  The  ratio  of  inner  aspect 
ratio  over  outer  aspect  ratio  becomes  larger  as  the  core  volume  frac¬ 
tion  becomes  smaller. 

Altogether  seven  variables  are  uBed  to  characterize  a  confocally 
coated  spheroid.  These  variables  are  inner  and  outer  complex  re¬ 
fractive  indices,  core  volume  fraction,  inner  aspect  ratio  and  outer 
aspect  ratio.  As  a  result  of  the  confocal  constraint  only  two  of  the 
last  three  variables  are  independent  leaving  Bix  independent  variables 
governing  extinction.  Projecting  contours  onto  a  two  dimensional 
plane  was  an  effective  technique  to  locate  strong  extinction  pro¬ 
duced  by  a  solid  spheroid  where  there  were  only  three  governing  in¬ 
dependent  variables  but  not  when  there  are  bIx  variables.  The 
appropriate  procedure  for  seeking  out  extinction  maxima  in  six  dimen¬ 
sions  is  to  employ  a  function  extrema  search  algorithm.  Such  a  search 
is  being  undertaken  and  it  is  too  early  to  describe  the  results,  how¬ 
ever  it  la  valuable  to  portray  in  two  dimensions  what  happens  to  the 
strong  extinction  obaerved  earlier  due  to  high  aspect  ratio  Bolid 
metal  spheroids  when  a  dielectric  confocal  coating  is  applied  over  the 
outaide. 

Calculations  have  been  made  for  metal  core  dielectric  coated 
spheroids  and  It  waa  found  that  dielectric  coatings  reduce  extinction 
at  all  aspect  ratios  not  close  to  one  for  both  oblate  and  prolate 
spheroids.  On  the  other  hand  it  has  been  discovered  that,  applying  a 
metal  coating  over  a  dielectric  core  can  significantly  increase 
extinction  for  both  prolate  and  oblate  spheroids.  Figures  6  and  7  map 
the  extinction  difference  between  coated  and  uncoated  oblate  apherbids. 
Extinction  due  to  solid  oblate  spheroids  with  real  index  20  and  imagi¬ 
nary  index  20  is  subtracted  from  the  extinction  due  to  coated  oblate 
spheroids  with  dielectric  core  real  index  1.5  imaginary  index  0.U1 
and  metal  coating  complex  refractive  index  equal  to  that  of  the  solid 
spheroid.  In  figure  6,  the  solid  spheroid  aspect  ratio  is  set  equal 
to  the  inner  aspect  ratio  while  in  figure  7  it  is  equal  to  outer  aspect 
ratio.  TheBe  contour  maps  are  similar  to  those  of  a  prolate  sphteroid 
except  that  the  contour  marking  the  boundary  where  extinction  of 
both  coated  and  uncoated  become  equal  occurs  at  an  aspect  ratio 
nearly  an  order  of  magnitude  smaller.  A  substantial  Increase  in 
extinction  i®  evident  in  figures  6  and  7  as  a  result  of  the  dielectric 
core  up  to  aspect  ratios  of  nearly  one  thousand  for  core  volume  frac¬ 
tions  between  0.1  and  0.9.  At  larger  aspect  ratios  the  solid  spheroid 
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•  is  superior,  however  this  cutoff  depends  on  the  optical  constants  of  ' 
the  metal. 

In  figure  8,  we  can  see  what  happens  to  this  coated  spheroid  at  an 
outer  aspect  ratio  of  one  hundred  and  core  volume  fraction  of  0.5 
when  we  adjust  the  core  refractive  index  to  determined  what  improve¬ 
ments  are  possible  over  the  dielectric  core  just  discussed.  The  orig¬ 
inal  dielectric  core  with  real  index  1.5  and  imaginary  index  0.01  is 
in  an  area  of  relatively  strong  extinction  but  figure  8  shows  that  by 
increasing  the  core  imaginary  index  above  0.01  extinction  will  get 
larger.  The  optimal  core  would  have  a  small  real  index  and  an  imag¬ 
inary  index  between  5  and  30;  values  to  be  expected  in  a  metal  a 
frequencies  just  below  the  plasma  frequency  at  visible  wavelengths. 

As  mentioned  earlier  the  Rayleigh  ellipsoidal  theory  is  not  a 
good  approximation  for  particles  that  hhVe  edges  and  vertices  such  as 
polyhedra.  Fortunately  dielectric  continuum  and  lattice  dynamical 
calculations  predict  the  absorption  properties  of  cubes  and  other 
polyhedra  small  compared  to  wavelength.  The  resiiltB  of  such  calcul¬ 
ations  for  a  cube  appear  in  figure  9.  Once  again  resonances  emerge 
along  the  K  axis  but  now  there  are  six  resonances  located  between 
imaginary  Indices  of  one  half  and  two  with  the  strongest  resonance 
appearing  furthest  up  the  axis.  Extinction  far  removed  from  these 
resonances  or  outside  the  restrahl  region  is  almost  indistinguishable 
from  that  of  a  small  sphere. 


Figure  9.  Extinction  isopleths  fo.  a  cube  as  a  function  of  complex 
refractive  index. 
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CONCLUSIONS 


When  particle  size  becomes  small  compared  to  wavelength,  extinc¬ 
tion  not  only  depends  on  geometric  cross  sections  but  also  becomes 
a  very  strong  function  of  complex  refractive  index.  Regions  of 
resonant  extinction  were  located  in  the  complex  refractive  index 
plane  for  a  variety  of  spheroidal  and  confoeally  coated  spheroidal 
particles  using  the  Rayleigh  ellipsoidal  approximation.  The  strongest 
most  extensive  is  found  for  high  aspect  ratio  metal  prolate  and  ob¬ 
late  spheroids  approximating  thin  needles  and  disks.  Even  stronger 
extinction  was  found  to  occur  when  metal  needles  and  dlBks  have  a 
dielectric  core.  On  the  other  hand  a  dielectric  coating  reduces  ex¬ 
tinction.  The  unbounded  peaks  of  extinction  resonances,  located 
where  the  real  part  of  the  complex  refractive  index  becomes  zero, 
were  proven  to  be  inaccessible  to  any  material  obeying  the  Kramers- 
Kroenig  relation  which  requires  absorption  at  restrahl.  Finally  the 
six  resonances  of  a  cubic  particle  were  explored  in  the  complex  re¬ 
fractive  index  plane.  Because  a  cube  ie  an  example  of  where  the 
Rayleigh  ellipsoidal  approximation  fails  most  severely,  the  extinc¬ 
tion  was  prndicted  using  a  result  taken  from  the  dielectric  con¬ 
tinuum  and  thb  lattice  dynamical  theories.^ 
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